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If the weak phase of B'^-B^ mixing (2/3s) is found to be significantly different from zero, this is a 
clear signal of new physics (NP). However, if such a signal is found, we would like an unambiguous 
determination of 2/3s in order to ascertain which NP models could be responsible. In addition, in 
the presence of NP, the width difference AFs between the two Bs mass eigenstates can be positive 
or negative, and ideally this sign ambiguity should be resolved experimentally. Finally, in order to 
see if the NP is contributing to in addition to M12, the precise measurement of |ri2| is crucial. 
In this paper, we consider several different methods of measuring B^-B^ mixing using two- and 
04 ' three-body decays with h cus and b ucs transitions. We find that the most promising of these 

is a time-dependent Dalitz-plot analysis of B^{B'^) D^pKK. With these decays, all of the above 
. issues can be addressed, and the measurement of the weak phase 7 is also possible. We also note 

that, with all three-body decays it is possible to resolve the sign ambiguity of AFs even without 
I determining CP phase (f)s. 

I. INTRODUCTION 

' Over the past several years, a number of discrepancies with the predictions of the standard model (SM) have been 
observed in B decays, intriguingly all in 6 — ^ s transitions. Some examples are: (i) in — >■ nK decays, it is difficult 
I to account for all the experimental measurements within the SM [l|, (ii) the values of the B'^-B^ mixing phase 

D , sin 2/3 obtained from different penguin-dominated & — >■ s channels tend to be systematically smaller than that obtained 
• from B^ — )■ J/ipKs 3, (iii) the fractions of transversely- and longitudinally-polarized decays in _B — !• (f)K* (/^ and 
] Jl, respectively) are observed to be roughly equal Q, in contrast to the naive expectation that /t/Jl ^ 1, (iv) the 
■ differential forward-backward asymmetry of leptons in the exclusive decay B — >■ K* fj,^ fj,~ is found to differ from the 

I> " SM expectations in both the low- and high-g^ regions {q^ is the dilepton invariant mass) @ . 

On . In light of this, it is particularly important to study 5 — ^ s transitions and look for new-physics (NP) effects. Now, 
if NP is present in AB = 1 b ^ s decays, it would be highly unnatural for it not to also affect the AB = 2 transition, 
1^ . in particular B^-B^ mixing. In order to see where NP can enter, we briefly review the mixing. In the Bs system, the 
I ■ mass eigenstates Bl and Bh {L and H indicate the light and heavy states, respectively) are admixtures of the flavor 
00 ] eigenstates i?° and B^: 
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\Bl) = p\B^)+q\m) , 

\Bh) = p\B^s)~<l\B's) , (1) 
with -I- |qp — 1. As a result, the initial flavor eigenstates oscillate into one another according to the Schrodinger 
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where M = Aft and F = F''' correspond respectively to the dispersive and absorptive parts of the mass matrix. The 
off-diagonal elements, — -Af|j* and Ffj = F|^, are generated by B^-B^ mixing. We define 

F, = ^"^^"^ , AM, = Mh " Ml, AF, = F^ - Fh . (3) 
Expanding the mass eigenstates, we find, to a very good approximation Q, 

AM, = 2|Mi^2l . 
AF, = 2|F^2|cos(/), , 

n 1 

(4) 
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where (j)s = aig{—M 12/^12) is the CP phase in Ai? = 2 transitions, and 2/3^ = arg(Mf2)- In Eq. (jj]) the small 
expansion parameter a is given by 
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12 



sin 0s . (5) 



This is expected to be ^ 1, and hence can be neglected in the definition of q/p. It is also important to note that the 
sign of AFs is equal to the sign of cost/)^, and in the case where there is no NP in T\2, the CP phase 0s = — 2/3s. 

The precise measurement of AMg determines \Ml2\ [!]• However, because of hadronic uncertainties, the SM 
prediction for AAfg is not very precise - in Ref. [9|, it is noted that the theoretical uncertainties still allow new- 
physics contributions to jMfjl of order 20%. In addition, Pfj can be calculated from the absorptive part of the B^-B^ 
mixing box diagram, leading to AFs. Unlike the Bd system, where AP^ is negligibly small, in the Bg system AP^ is 
expected to be reasonably large, which leads to certain advantages for the search for CP- violating effects in the Bs 
system over that of Bd system. The updated SM predictions of the width difference and the CP phase 0^ are given 

by i 

APf*^ ~ 2|P^2l = 0-087 ±0.021 ps-\ 

0s ~ 0.22° . (6) 

Although the SM predictions for AAfs and APs are not precise, the SM does predict that 2/?s — 0, which makes 
it a good observable to use in the search for NP. Consider first the case where the NP contributes only to Ml2- If 
one measures a value for 2/3s that is significantly different from zero, this will indicate NP in B'^-B'^ mixing (in Ml2)- 
However, to cover all bases, one more step must be done. Suppose that NP is present, but it produces 2/3s — 180°. 
Now indirect CP violation, which measures sin2^s, will not give a signal. But one can still detect the NP by measuring 
the sign of APs - in the SM, APs > 0, while it is < if 2/Js = 180°. Also, even if NP is discovered through indirect 
CP violation (i.e. sin2/3s 7^ 0), this only determines 2/3s up to a twofold discrete ambiguity. Since the sign of cos2/3s 
can be determined by the sign of APs , the knowledge of this sign is one possibility to remove this discrete ambiguity. 
Alternatively, one may try to find a method which allows a direct determination of 2/3s without any ambiguity. Now 
suppose that the NP also contributes to Pf2. Since there is NP in A/fj, its presence can be detected as above. But now 
the twofold ambiguity in 2/3 s cannot be removed from the knowledge of the sign of APj, since this only determines 
the sign of cos 0s (and not cos 2(3s). Thus, for the case where there is NP in Pf2, one must find another way to remove 
the twofold ambiguity in 2/3s. 

The CDF and D0 11 1 collaborations have measured the CP asymmetry in B'^ J/'4'4>, and found a hint 



for indirect CP violation. In general, this result is interpreted as evidence for a nonzero value of 2pf't', and the 
contributions of various NP models to the Bs mixing phase have been explored [T2| - [l8j . It has also been pointed 
out that NP in the decay b — sec could also play an important role [l^. Recently CDF and D0 updated their 
measurements of the CP-violating phase. The 68% C.L. allowed ranges are [20I [2l| 

2^s^"^ e [2.3°, 59.6°] U [123.8°, 177.6°] , CDF, 

€ [9.7°, 52.1°] U [127.9°, 170.3°] , D0 . (7) 

Most of the values of 2/3^"^ here suggest NP. 2/3^"^ is obtained with the twofold ambiguity 2/3^"^ o tt - 2/3^'*', and at 
present there is no preference for either of the two solutions. As mentioned above, the possibility of NP in the decay 
h sec cannot be ruled out, so that the phase 2/3^*^ extracted from B^ — > J/'00 should not necessarily be taken as 
purely a mixing phase. It is therefore worthwhile to look for a process in which NP in the decay can essentially be 
neglected, and which permits the determination of 2/3s without any ambiguity. If the measured value of 2/3s is found 
to be significantly different from that in B^ J it will be clear signal of NP in 6 — >■ sec. 

In addition, the D0 Collaboration recently found a large CP asymmetry in the like-sign dimuon signal, which they 
attribute primarily to a|.^, the semileptonic CP asymmetry in B^ — Xg^iv [13, H^. Now, the D0 result is less 
than 2a away from zero and consequently to an excellent approximation also about 2cr away from the SM prediction 
(^-SL*^ ~ 2 X 10~^) [9]. Still, NP in B^-B^ mixing can explain the result (for exarnple, see Ref. J4]). However, if one 
wishes to reproduce the central value of a^^, one requires NP specifically in Ff2 [2^, [2^. There are NP models that 
can contribute to F^2 through the decay b — >■ sr+r" [13, H^, and a significant enhancement of its magnitude over 
that of the SM M is possible. Furthermore, the possibility of NP effects in TI2 through the decay b sec cannot be 
ruled out [H ■ 

We therefore see that there are some hints of NP in the Bg system, but nothing definitive yet. Thus, it is important to 
look for additional methods of probing NP in B'^-B'^ mixing. Ideally, the new method(s) would allow an unambiguous 
determination of the mixing phase 2/3s. Also useful are methods which remove the sign ambiguity in APs even without 
providing any direct information on the CP phases 0s or 2/3s. Finally, if NP is present in the mixing, we would like 
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to know if it contributes to Ffj in addition to Mfj- Hence, along with the removal of the sign ambiguity in AF^, 
independent and unbiased measurements of iTfjl and are essential. 

Several years ago, the two-body decays — t- DfK^,D*^K^,... were examined with the idea of extracting 
weak phases [1^. Because the final state is accessible to both and mesons, a mixing- induced indirect CP 
asymmetry occurs. Using this, and assuming that AF^ is sizeable, the conclusion of Ref. [29] is that one can measure 
the phase 2/3s + 7 with a twofold discrete ambiguity, and that this ambiguity can be removed if factorization is 
assumed. However, if there is NP in mixing, AF^ < is allowed as well. This implies that, in fact, 2/3s -I- 7 

can be obtained with a fourfold discrete ambiguity (or twofold if factorization is assumed). 

In Ref. [30] it was shown that the sign ambiguity in AF^ can be removed using B° ->■ DfK"^ decays. Although 
the method does not allow a direct determination of the phase 2/3<,, it does discriminate between the two solutions 
with cos2/35 > and cos2/3s < 0, which then determines the sign of AFg. However, the method is based on several 
assumptions: (i) the weak phase 7 is taken from the _B-factory measurements, (ii) factorization is assumed, i.e. the 
strong phase is taken to be ~ 0, and (iii) the SM-predicted value of Ffj has been used in the analysis. 

In 1991, the decays B^{B^) — >• D^pcj), where D^p is a neutral D-meson CP-eigenstate, were proposed to extract the 
CKM angle 7 with a twofold ambiguity [Sllls^l- However, these methods assumed that the phase 2/3^ is approximately 
zero (or known) . The current experimental data [see Eq. ([7])] is not completely in favor of this assumption ~ there is 
the possibility that 2/3^ can be significantly different from zero. In addition, at present 2f3s is measured with a twofold 
ambiguity, which adds a further discrete ambiguity to the determination of 7. 

We therefore see that previous analyses of two-body B decays only partially probe NP in B^-B^^ mixing ~ 2/3s is, in 
general, not determined unambiguously, the sign ambiguity in AFg is generally unresolved, and the possibility of NP 
affecting has not been considered. In this paper we go beyond the previous analyses to explore all of these issues. 

In Sec. in] we review the two-body decays. In particular, in Sec. Ill C[ we update the analysis of B^{Bl) — )• D'^pcj), 
considering both AF., > and AFs < 0. In Sec. IIIIl we present the Dalitz-plot analyses of three-body decays. In 
particular, in Sec. IIII C[ we focus on B",{B°) OfK^-K^^Df-K^K^,..., using the interference between the different 
intermediate resonant decays to provide additional information. And in Sec. IIIIDi we show that a much greater 
improvement can be obtained by performing a time-dependent Dalitz-plot analysis of the decay (5°) — >■ D^pKK. 
Finally, in Sec. lIVl we present a possible way to determine AF^, or equivalently jFfj] and (ps, using three-body decays. 
We conclude in Sec. |Vl 



II. TWO-BODY DECAYS 
A. B°iB°)^f,f 



Consider a final state /, not necessarily a CP eigenstate, to which both and can decay. In the presence of 
B^-B^ mixing, the time-dependent decay rates are given by [33j 



r(i?°(t)->/) 



r(s°(t)->/) 



ie-r=*{(lA/l2 + \Aj\^)coshiAT,t/2) + {\Af\^ - \Aj\^) cos Am^t 



2sinh(AFsi/2)Re 



2sinAmoiIm 



} 



ie-r=*{(lA/l2 + \Af\')cosHArst/2) - {\Af\' - \Aj\') cos Am^t 



2sinh(AFst/2)Re 



2 sin Amoi Im 



'-A}Af 



} 



where Af = A{B° -y f), Af = A{B°^ -> /), and q/p = e-^*^^. This yields 



F(B,"(t) ^ /) -F(B,"(t) ^ /) - {\Af\'^ + \Af\'^)e-^'*[C cos Arust- S sin Arust] , 
F(B°(t) ^ /) +F(SO(t) ^ /) ^ (JA/l^ + lA/l2)e-r=* [cosh(AF,V2) - -4Ar sinh(AF,,i/2)] 



where 



1_^ 2ImA _ 2ReA 

" I + IAI2 ' I + IAI2 ' I-HAI2 



pAf 



(8) 



(9) 



(10) 



The idea is that, by fitting the data corresponding to the difference ("tagged") and sum ("untagged") of decay rates 
to the four time-dependent functions given on the right-hand side of the equations in Eq. (jH]), the coefficients of these 
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functions can be obtained, from which C, S and Aat can be derived. However, there is a comphcation - in the 
presence of NP in AB — 2 transitions, AF^ is unknown (though it is assumed to be reasonably large). Therefore, for 
the untagged combination, both ATg and Aat must be found in the fit. Still, though this will determine lAFs], its 
sign will remain unknown. The reason is that only the function smh{ATst/2) is sensitive to the sign of AFs, and it is 
multiplied by Aaf- Thus, any change in the sign of AFs can be compensated for by changing the sign of Aaf- The 
bottom line is that any analysis which uses Aaf will have a discrete ambiguity due to the unknown sign of AFs. 
Similarly, 



nm) ^ I) 



1 



2sinh(AFst/2)Re 
2sinh(AFst/2)Re 



Aj\^) cosh(AFsi/2) + (l^jp - | A/p) cos Am^t 



- A*-A- 



Aj\^)cosh{ATstl2 



2sinAmsiIm 
(1^/ 



2 I J J2 



- A*-A - 



2sin Amcilm 



AjAf 



where Af = A{B"^ f) and Aj = A{B° /). Then 

T{B"At)^f)-rmt)~^f) 



C cos Arrigt + S sin Amst 



where 



T{BO{t) f) + T{B°{t) ^ f) cosh(AF,V2) - ^Ar sinh(AFsi/2) ' 
|A|2 ^ 2ImA 



C = 



1 



1 + |A|^ 



S = 



1 + IAP 



-4Ar = 



2ReA 



qAf ■ 



(11) 



(12) 



(13) 



B. B%m)^DfK^ 

Consider the decay B1{B1) — > PP {P is a pseudoscalar), in which the final state contains a single c quark^. 
Excluding those final states involving jy's, there are only two decays in which the B^ and B^ amplitudes are of 
comparable size: 

B°s{B°s)r^ D7K+ and S^(B°) D+K'. The decays are mediated by color-allowed tree-level 
transitions b — cus and h — > ucs. Within the SM, the amplitudes take the form^ (there is a minus sign associated 
with the u quark) 

A(B," ^ D-K+) = r, A(i?o ^ DtK-) = -f , 

A{B° D-K+) = f'e-'^ , A(B," ^ DtR-) = -T' . (14) 

We have explicitly written the weak-phase dependence, while the diagrams contain strong phases. The magnitudes of 
the Cabbibo-Kobayashi-Maskawa (CKM) matrix elements iV^^KisI and |Vu*bV"cs| have been absorbed into the diagrams 
T' and T', respectively. (As this is a & — > s transition, the diagrams are written with primes.) 
Using the amplitudes of Eq. (fT4| . one obtains [see Eqs. (fTO|) and (flB])] 



^=TW' S = -j^sH^Ps+^-5)^ ^^, = _|^cos(2/3,+7-5) 



^ " TT^ ' ^"T+W™^^'^''^^^'^^ ' -^^r = y|^cos(2/?.+7-|-<5) , (15) 

where |A| — T' /T' (defined to be positive) and 5 is the strong-phase difference between T' and T' . |A| can be obtained 
from the measurement of C. Using this, S and Aav give sin(2/3s -1-7 — 5) and cos(2/3s -f- 7 — 5), respectively. Thus, 



^ Much of the discussion in this subsection can be found in Ref. :2£|, except that here NP in AFs is considered. 

^ In Ref. |26|| . it is shown that NP in the decays b — > cus and b — > ucs is strongly constrained. Such NP contributions are therefore 
neglected throughout this paper. 
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one obtains 2/3^ + "f — S with no discrete ambiguity. Similarly, 2/3^ +7 + 5 can be obtained with no discrete ambiguity 
from S and Aaf- These can be combined to give the phases {2(3g + 7, 6) with a twofold ambiguity [(2/3^ + 7, S) or 
(2/3s + 7 + TT, (5 + tt)]. This discrete ambiguity can be removed if one assumes factorization, which predicts S to be 
near 0. 

In fact, this is not quite correct. As discussed below Eq. (flOl) . in the presence of NP in B^-B'^ mixing there is an 
additional discrete ambiguity due to the unknown sign of AFs. Thus, the two-body decays B^{B^) — >■ DfK^ permit 
the extraction of 2/3^ + 7 with a fourfold ambiguity (or twofold if factorization is assumed). 

Now, the value of 7 can be taken from the independent measurements at the i3-factories. One then obtains 2/3s 
with a fourfold ambiguity. Alternatively, since 7 has not been measured in decays, it can be kept with the aim 
of determining its value independently (this was the original purpose of Ref. [23] ■) We adopt this latter approach in 
much of the paper. 

We therefore see that this method permits the extraction of 2/3s + 7 with a fourfold ambiguity (or twofold if 
factorization is assumed). It does not resolve the sign ambiguity in ATg, and says nothing about the possibility of 
NP affecting Ffj- In order to address these remaining points, it is necessary to examine other methods. A first step 
involves the decays 5° (i?^) — > D(jp(j>, discussed in the next subsection. 



C. B°iB°) D^p4> 

Another pair of decays to which the method of the previous subsection can be applied is (i3°) D^4>- Here 

the decays are mediated by color-suppressed tree-level transitions. The amplitudes (of comparable size) are given by 

A(i?° ^ i?V) = -Cfe^^ , AiB'l -> 5» = Ct , 

A{B° ^ ^» = Cf e^^T , A{B"g ^ D°^) = -C^ ■ (16) 

By measuring the time dependence of the decays, one can obtain S*, Aaf and Aaf as given in Eqs. (ITU|) and (|13l) . 
Using these observables we define 

sin(2/3, +-f + 5^) = -^^1^^ = Sd , sin(2/3, + - S^) = = So , 

cos(2/3, +j + 6^) = ili^^Ar = r , cos(2/?, +j~6^)^ i±^i^P = if r , (17) 

with — arg(Cf/C2 ). The method of the previous subsection then allows us to obtain 2/3^ + 7 with a twofold 
ambiguity (for the moment, we put aside the ambiguity due to the sign of AFs). 

The advantage of these decays is that there is a third decay which is related: i?°(i?°) — )• D^p(f), where D^p is a CP 
eigenstate (either CP-odd or CP-even). In our analysis we consider D^p as the CP-even superposition {D^ -\- D'^) / \/2. 
The amplitudes for the decays are then given by 

V2A(B," ^ p0) = -Cf e'^ + ct , 

V2A{b1 D'ijpc^) = Cfe-'^' - ct ■ (18) 

By measuring the time-dependent decay amplitudes of B^{B'^) -> Dcf) {D = D'^ , D'^ , D'^p), one can extract the 

magnitudes \Cf\, |C||, \Adcp \ = \A{B^ D^p<P)\ and \Adcp \ = \AiB° D^p(t>)\ (they are combinations of the 
overall normalizations and the C parameters [Eq. ((TU)) ]'). 
Using the first equation of Eq. (flSj) , we define 



^ 2\Ct\\Ct\ 



Similarly, from the second equation of Eq. (jlSp . we get 



* 2\Cf\\Ct\ 

Therefore, in the case of the B'^{B'^) D(j) decays, we have two more observables, E+ and E~. Combining Eqs. (|17l) . 
(IT^ and ([^0]) . it is straightforward to find expressions for sin2/3s, cos2/3s, sin(2/3s + 27) and cos(2/3s -I- 27) in terms 
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of the above observables: 



sin 2/3s 




2 



sin(2/3, + 27) 




2 



cos 2/3s 




2 



cos(2/3s + 27) 




(21) 



with 




sin2/3s and sin(2/3s + 27) could be obtained, which correspond to determining 2/3^ with a twofold ambiguity and 27 
with a fourfold ambiguity. In the present case, the dependence on AF^ is included. This allows us to obtain A^p and 
A£p, which then permits us to measure cos 2/3^ and cos(2/3s + 27), in addition to sin 2/3^ and sin(2/3s + 27) [Eq. (^1]) ] . 
These measurements allow an unambiguous determination of 2/3s and 27. We therefore see that a nonzero AT^ helps 
quite a bit in determining the weak phases. As has been discussed above, the sign of AFg is not known, which implies 
that and also have a sign ambiguity. This means that, in fact, 2/3s and 7 are determined up to a twofold 
and fourfold'^ ambiguity, respectively. Therefore, once we are able to fix the sign of AFs, the Bg{Bg) -> Del) decays 
might be considered as an alternative mode to probe simultaneously 7 and 2/3s. 

We therefore see that two-body b — cus/h — > ucs decays do not provide sufficient information to measure the CP 
phases 2(3 s and 27 in an unambiguous manner. In the next section we show that there are several ways to improve 
upon the two-body decay methods by using a Dalitz-plot analysis of the corresponding three-body decays. 



In recent years, it has been shown that one can get useful information from three-body B decays. For instance, 
time-integrated Dalitz-plot analyses of B^ — )■ K-ktt and 5° nKK decays have been proposed as a probe of 7 (s^ . 
And various tests of the SM, as well as the extraction of weak phases, have been examined in the context of B — )■ Kinr, 
B KKK, B nKK and B tttttt decays [||. 

In the previous section we discussed two-body b — >■ cus/b — > ucs decays; in this section we examine the corresponding 
three-body decays. In _B^(i?°) — >■ PPP decays which receive a tree contribution, there are 5 final-state (/, /) pairs: 



{D^K+y, D+K-TT°), {DjK°TT+, D+K°TT-), {D-K+K°, D+K°K-), {D°K+K', D^K+K~), and {D°K°K°, 



D'^ K'^ K^) . The CKM matrix elements of these decays are the same as in the corresponding two-body decay modes, 
and will therefore exhibit very similar time-dependent CP asymmetries. 

The decay amplitude of 5^(5") — ?► PPP receives several different contributions, both resonant and non-resonant. 
In the following, we perform a time-dependent Dalitz-plot analysis of the three-body decays. This permits the 
measurement of each of the contributing amplitudes, as well as their relative phases. As we will see below, the 
Dalitz-plot analysis reduces the ambiguity in the measurement of 7 and 2/3^ compared to the corresponding two-body 
decays. We also show how this analysis resolves the sign ambiguity in AF^. 



Here we review certain aspects of the Dalitz-plot analysis. We focus on the general three-body decay B — > PiP2^3- 
We define the Dalitz-plot variables 



^ Using Eq. II21I I. we can determine cos 27 without any ambiguity, whereas, due to the unknown sign of j4^p or ^^p, sin 27 can be 
determined only with a twofold ambiguity. Combining these two results, 27 can therefore be determined with a twofold ambiguity (or 
7 with a fourfold ambiguity). 



III. THREE-BODY DECAYS 



A. B°{B°)^f,f 



B. Dalitz-plot analysis 



S12 = 



bl+P2)^, Si3 = (pi +P3)^ , S23 ^ {P2 + Ps)"^ , 



(23) 
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which are related by the conservation law 



Sl2 + Sl3 + S23 = mi 



- rrin 



(24) 



This shows that there are only two independent variables (below, we use si2 and S13). 

B — > P1P2P3 can take place either via intermediate resonances or non-resonant contributions. A widely-used 
approximation in the parametrization of the decay amplitude is the isobar model. In this model, the individual terms 
are interpreted as complex production amplitudes for two-body resonances, and one also includes a term describing 
the non-resonant component. The amplitude is then written as 



-^(S12,S13) = y^^ajFj{si2,si3) 



(25) 



where the sum is over all decay modes (resonant and non-resonant). Here, the aj are the complex coefficients 
describing the magnitudes and phases of different decay channels, while the Fj[si2, S13) contain the strong dynamics. 
The CP-conjugate amplitude is given by 



^(S12,S13) = y^ajfj(si3,Si2) 



(26) 



where Fj{si3,si2) = Fj{si2,si3). 

Now, in the experimental analysis, the Fj(si2, S13) take different (known) forms for the various contributions. By 
performing a maximum likelihood fit over the entire Dalitz plot, one can obtain the magnitudes and relative phases 
of the Uj, and similarly for the dj. Thus, the full decay amplitudes can be obtained. 



C. B"^{B°)^D^K^^° 

In this subsection we focus specifically on the decay B'^{B'^) — > DfK^n^, and use a modification of the method 
elaborated in Ref. [3()|. The Dalitz-plot variables are 



The amplitudes are written as 

^(s+,s-) =^aj.Fj(s+,s-) , Ais+,s-) = Y,ajF,{s-,s+) 



(27) 



(28) 



The time-dependent decay rates for the oscillating B^{t) and B'^{t) mesons, decaying to the same final state /, are 
given by 



ns^sit) ^ /) 
nmt) ^ f) 



Ach(s+, s") cosh(Ar,t/2) - A,h{s+,s-) sinh(Ar,t/2) 
-I- ylc(s^, s^) cos(Amsi) — yls(s+, s^) sin(Amst) , 



Aphis', s+) cosh(Ar,t/2) - Ash{. 



)sinh(Ar,t/2) 



Ac{s , s+) cos(Amst) + As(s , s^) sin(ATOst) 



Here 



(29) 



Ach{s+,S-) = |^(S+,S-)|2 + |^(S+,S-)|2 , 

A,{s+,s-) = \A{s+,s-)\^-\A{s+,s')\^ , 
Ash{s+.s-) = 2-Re{e-^'P^A{s+,s')A*{s+,s~)) , 
A,{s+,s-) = 2Im(e-2'''=^(s+,s-)yl*(s+,s-)) , (30) 

where ^(s+, s~) and ^(s+, s~) are given in Eq. (^5)) . 

Now, there are a number of resonances which contribute to B^{B^) DfK^n'^. For illustrative purposes, we 
consider just two of them: DfK*^ (892) and Dl^K^. The decays K*^ K^tt^ and D;^ ^ Dfir^ have already 
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been_ observed: B{K*^{892) K^tt°) = 50%, ^ Dfir") = (5.8 ± 0.7)% 50(5°) ^ Li±7^*=F and 

BO(BO) ^ D^^ifT are the 50(^0) ^ equivalents of the decay discussed in Sec. IHl B°{B°) DfR"^. The 
additional ingredient here is that we also consider the decay products of the V, so that we have the full decay chain 
50(5°) ^ VP^PPP. 

For these two resonances, we have 

Ak'[B^s 
Ad'AB", 



D-K*+ ^ D;K+7t") = afe'-'FK^ , 
D-K*+ D;K+Tr") = afFK' , 

D*-K+ ^ D-K+tt") - 
Including both resonances, the amplitudes of 5^(5^) — ?> Dj K^n^ are 



a^' I'D' 



A = e'^afFK.+af'Fo.) , 
A = (af>A- +af=>D*) . 
With these amplitudes, Ach, Ac, Ash and Ag [Eq. (1501) ] take the forms 



^ch 



(i«r 



K' |2 



FK-A''+{\aT^? + \aT^?)\Fot? 
2Re ((af*FK.)*(af'-F^i?;)) + 2Re (^{af Fk-)* {af Fd* 
(l«r 1^ - l«f 1^) \Fk^ ? + (|af \' - |af 1^) IP,,, P 



+ 2Re ((af FK.)*(af''^i3:)) - 2Re (^{af Fk^)* {a^^ Fd-^ 
cos(2/3, + 7 + <5k. )|af' I |af ' I I^K- P + cos(2/3, + 7 + (5i5. )|af '* | |a, 
+ cos(2/3, + 7 + 5)|af' ||af |Re {F^,Fd'J 
- sin(2/3, + 7 + (5)|af * ||af = |Im {F^,Fd'J 



Rc 



-sin(2^^ + 7 + (5K.)|af 



K* \\„K' 



\Fk'\''- sin(2/3, + 7 + (5i3.)|af ° 1 1": 



Fr 



- sin(2/3, +7 + ^)|af ||a^= |Re (P^.Po* j 
+ cos(2/3, + 7 + 5)|af * ||af |Im {F^.Fd-J 
Im e-(2^=+^)(af Fz,.)*(af>KO)' 



(31) 



(32) 



(33) 



with Sk* = - arg ((af *)*af *), ^d* = - arg (^(af= )*af"^, and (5 = - arg |^(af *)*af 

Above, in the discussion of the time-independent Dalitz-plot analysis, we noted that the magnitudes and relative 
phases of the aj can be obtained from a maximum likelihood fit over the entire Dalitz plot, given assumed forms for 
the Fj^s. The same holds true for the time-dependent Dalitz-plot analysis - the magnitudes and relative phases of the 

contributing resonances, i.e. af^ , af^ , af ° , and , can all be obtained. Indeed, such an analysis has been performed 
by the Babar and Belle collaborations for the decay B^{t) — ^ Ks'n'^n~ [s^]. In particular, all the coefficients that 
multiply the F*Fj [Eq. p3p ] bilinears can be obtained from a maximum likelihood fit to the corresponding Dalitz 
plot PDFs. 

This permits the extraction of the weak phases. For example, we can extract 2/3^ -1-7-1-5 without any ambiguity 
from the third and fourth terms of Af"^"^. In a similar manner, the time-dependent Dalitz-plot analysis of B'^{B'^) — > 
DfK~'!T° allows the extraction of the phase 2(3s + j — S. The combination of these two results yields 2/3^ -I- 7 and 
S with a twofold ambiguity. And if factorization is imposed, the discrete ambiguity is removed entirely (only the 
solution with 5 ~ is kept). The key point here is that we do not use A^/^ ^ at all. As a consequence, there is no 
discrete ambiguity due to the sign ambiguity of AF^ [see the discussion following Eq. ([50)) ]. This is to be contrasted 
with two-body decays. There 2/3^ -I- 7 can also be obtained with a twofold ambiguity. However, because ,4Ar and 
^Ar are used [Eq. ([T5|)]. there is an additional discrete ambiguity due to the unknown sign of AF^. 
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We note that one can extract different trigonometric functions such as | sin(2/3s + 7 + | cos(2/3s + 7 + 
I cos(2/3s + 7 + di)\, etc., from A^f^^'^ [Eq. (l33l) ]. Due to the sign ambiguity of AFj, which can be viewed as the sign 
ambiguity in A^y^ ^ , the sign of these trigonometric functions cannot be determined. Depending on the sign of AF^ , 
their sign could be positive or negative. Therefore, we can determine the sign of AF^ if we are able to fix the sign 
of these trigonometric functions. Now, the functions sin(2/3s + j + S) and cos(2/3s + j + 6) can be extracted without 
ambiguity from A^=^^ , which fixes the sign of AF^ and hence removes the discrete ambiguity in Af^^^^ . Note that 
this can be done without measuring (/is . This method can therefore be used to determine the sign of cos (f)s . 

In the above, we have concentrated on the decay 5^(5^) DfK^n^. However, any of the decay pairs discussed 
in Sec. |lll] can be used. All that is necessary is that there be at least two resonances contributing to the decay. We 
therefore see that, by using such three-body decays, one can obtain 2/3s +7 (with a twofold ambiguity if factorization is 
not assumed), as well as resolve the sign ambiguity in AFg. The resolution of the AF^ sign ambiguity determines the 
sign of cos (j)s . The precise knowledge of 7 from other measurements allows one to obtain 2/3^ with a twofold ambiguity 
(since 2/3s + 7 can itself be extracted with a twofold ambiguity), which can be compared with the measurement of 
2P, from 5° ^ J/^jcj) [Eq. ©]. 

Still, it is preferable to have a method that allows the direct determination of 2/3s and 7 individually. This can be 
done by measuring the decay B'^{B'^) — > D^pKK, which is discussed in the next subsection. 



D. ^ D^cpKK 

In Sec. Ill CI we discussed the two-body decays B^[B^) — > Dcj) [D = D°,D°,Dqp), and showed that it is possible 
to extract 2/3s and 27 with a twofold ambiguity due to the unknown sign of AF^. The time-dependent Dalitz-plot 
analysis of i3"(i3") — > D^KK,D^KK is similar to that of the previous subsection, with the intermediate resonances 
0(1020) or /o(1500) decaying to the final state KK. In this subsection we consider in addition the related three-body 
decays _B°(i3°) ^ D°cpKK, with D^p = 1/V2(D" ± D^). 

Bs(i?°) -> DqpKK receives contributions from several different intermediate resonances: 0(1020), 0(1680), 
/o(1500), /o(1710), D*f, etc., which follow the decay chains B°{B°) D°cp(f> D^pK+K-, B°(B^) ^ D^pfo 
Dl;pK+K-, B°{B°) ^ DlfK"^ D^pif^if^. To simplify our analysis, we consider only the 0(1020) and /o(1500) 
resonances. The amplitude with an intermediate resonance can be written as 

V2A^(B° ^ D%p(j}{-^ K+K-)) = A(B^ ^ D°K+K-) + A(B," ^ D^K+R-) , 

y/2A^{B° ^ Dl.p(t){-^ K+K-)) = A(B^ D^K+K') + Ai^B'l D^K+K') , (34) 

where 

yl(B° ^ D"K+K-) = -Cfe'''F^ , 

A{B° ^ D"K+K-) = -C^F^ , 

A(B° ^ ^> ^ D°K+K-) = CtF^ , 

A(B° ^ L'°0 ^ D"K+K-) = Cf e-^^F^ . (35) 

The amplitude with an intermediate fo resonance is given by a similar expression, with the replacement (j) fo. 
Including the contributions from these two resonances, the total amplitude can be written as 

^(5° ^ D%pK+K-) = A^iB'l ^ D"cpK+K-) + Af,{B", -> D"cpK+K-) , 

^(5° ^ D%pK+K-) = A^(B," ^ DlpK+K-) + %(B0 ^ D^pK+K-) . (36) 

With these, A^,^'^ and can be computed similarly to Eq. First, we have 

Ac'''' = E [(l^^l' - l^^l') + 2Re(A^A;.„ - A^Al)] , 

^f/"" = E + \A,?) +2Ke{A^Al+ A^A}^)] , (37) 

i='t>Jo 
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in which 



Re{A^A%-A^A%) = \Cl''\\Ct\smj[r^{sm6^Re{F^Fl) + cosS^lm{F^Fl)} 

+ r/o { sin df„ ReiF^F}^ ) - cos 6f„ lm{F^F}^ ) } ] , 

Re{A^A% + A^A%) = \Ci"\\Ct\[ReiF^Fl) - cos-f{ cosS^Re{F^Fl) - sin5^Im(F^F;j} 

- cosj{ cosSf„Re{F^FfJ + sin 6 f^lui{F^Ffj} 
+ r4>rfo{cos{S^ - 5/jRe(F^F;j - sin((50 - 5f„)lm{F^Fl)} 



(38) 



where n = \Ci\/\C^\ and Si = arg(Cj/Ci) (i = 0, /o)- Using Eq. ^ in Eq. a maximum hkehhood fit to the 

Dahtz-plot PDFs ahows one to extract 



|C|°||C||r,cos7 cos5, = , 
ICi^WC^lnsin-f cos 5, = al . 



This gives the ratio 



-4- = tan 7 



(39) 



(40) 



Since the hadronic uncertainties cancel in the ratio, it yields a theoretically clean determination of the angle 7 with 
a twofold ambiguity, even without the knowledge of the strong phases. 
Second, we have 



= Im [e-'^^^A*A] = Im [e''^^^ {A;A^ + A;Af, + A%A^ + AlAf,)] 
The first and fourth terms of A^^^ are given by 

Im [e-^'^^A*A,] = hm [e-^'^^ |C^nF,|2(l + rfe-^^T + r,(e*(^*-^) + e-'(**+'^))) 



(41) 



(42) 



which allows the extraction of sin 2/3s, sm{2/3s + 2^), sm{2(3s+j — 6^/fg) and sin(2/3s +7 + (50/jp). The 0-/o interference 
terms are given by 

Iu,[e-'^^^{A;Af,+A}^A^)] = llm[e-'^^^\Ct\\Ct\{~{F;Ff,+FlF,) 



rjo {e-^^''-''o)F;Fj,+e-^^''+'fo)FlF^ 



.( 



g-»(27+'5*-'5/o) F^Ffg + e^*(2^^''*+''. 



(43) 



This yields 



lm[e~'^^^{AlAf,+A%A^)] = [sin2/?, {Re(F;F/J + Re(F;„F^)} 



- sin(2/3, + 7 + S^)Re{F;Ff,) - cos(2/3, + 7 + S^)lm{F;Ff,) 
+ sin(2/3, + 7 - ^^)Re(F;^F^) - cos(2/3, + 7 - 6^)lm{FlF^)} 

- Tf, { sin(2/3, + 7 + <5/jRe(^^;„F^) - cos(2/3, + 7 + <5/JIm(F;„F^) 
+ sin(2/3, + 7 - 5f,)Re{F;Ff,) ~ cos(2/3, + 7 - <5/JIm(F;F;„) } 

+ r/o{ sin(<50 - Sf, + 2/3, + 27)Re(F;F/„) ~ cos(<5^ - <5/„ + 2/3. + 27)Im(F;i^/„) 



+ sin(,5/, ^5^ + 2/3, + 27)Re(F;„F^) - cos(% - 5^ + 2/3, + 27)Im(F;^F^ 



(44) 
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From the above, we can extract 



r, 



sin(2/3, +j±5i) = SJf^^ , cos(2^, + j ± Si) = C'^kk ' sin2/3s = Sdkk 



r„- sin(2/3, + 27 ± <5,,) = S^^i, , -n, cos(2/3, + 27 ± <5,,) = C]^±^ , (45) 

where r^j = and the corresponding Sij = Si — Sj = 4>,fo)- It is straightforward to find expressions for 
tan(2/3s + 7) and tan(2/3s + 27) in terms of the above observables: 

tan(2/3. + 7) = -^^^^-^^ , tan(2/3, + 27) = + ■ (46) 

'■-'DKK ^ '-'DKK ^DKK^^DKK 

With these, one can obtain the expression for tan 7 in terms of the extracted observables: 

tan(2/3, + 27) - tan(2/3, + 7) 
tan7 — . (47 

' 1 - tan(2/3s + 7) tan(2/3, + 27) ^ ' 



This way of getting tan 7 uses Af^^ [see also Eq. (|40l) ]. 
Combining Eqs. P5|) and (HT)) . we obtain 

tan2/?.^ tan(2/3 +7)-tan7 ^ 

l-tan(2/3,, +7)tan7 ^ ^ 

This determines 2/3s with the twofold ambiguity 2/3s tt + 2/3^. However, as we note in Eq. (H5|) . we can extract 
sin2/3s without any sign ambiguity. This determines 2 /3s with the twofold ambiguity 2/3^ — > tt — 2/3^, which is different 
from that obtained in tan2/3s. Therefore, the combined measurements of ta.n2/3s and sm2(3s allow us to extract 2f3s 
without any ambiguity. The sign ambiguity in AF^ can be resolved in a similar way to that discussed in Sec. IIII CI 

Above, we discussed the interference between the two resonance states 0(1020) and /o(1500). However, the anal- 
ysis would hold equally for the interference between any two resonances decaying to the same final state. Similar 
information can also be obtained from the time-dependent Dalitz-plot analysis of 3^(13^) — >■ D^pK^K^ . 



IV. EXTRACTION OF 4^ AND \Ti-^\ 

In the previous section(s) we examined methods for extracting the CP phase 2/3s using various two- and three-body 
decays. The idea is that if a nonzero value of 2/3s is found, this will be clear evidence of new physics in B'^-B'^ mixing. 
In addition, if such a value of 2(3^ is obtained, we will want to know its exact value in order to ascertain which different 
models of NP could generate such mixing. To this end, the best method will be that for which the discrete ambiguity 
in 2f3s is minimized. However, there is one question which has not yet been addressed: if NP in the mixing is found, 
does it contribute to Ff2 in addition to M;f2? 

In order to answer this question, (f>s and |rf2| must be measured and their values compared with those predicted 
by the SM [Eq. As pointed out by several authors, the possibility of NP in F^2 cannot be ruled out, and NP 

can enhance rf2 significantly above the SM prediction, even including the error bars. In order to measure or extract 
both (ps and |F^2l; '^6 need (at least) two observables which are sensitive to NP effects in F^2- We choose AFg and 
the semileptonic asymmetry a^j . As we show below, the precise measurement of these two observables can be used to 
extract 0s and |Ff2|- 

The expression for AF^ is given in Eq. the semileptonic asymmetry is defined as 



Im 



^ 12 



2^ 
AM, 



(49) 



Combining Eqs. (HJ and p9|) we obtain 



tan (p, = 



IF?,! = 



<i AM, 
AF, 



AF, 



,s 2 



AM, 



(50) 
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Now, AMs is known very precisely - AMs = 17.77 ± 0.12 [3,1331 - and so the precise measurements of a^; and ATg 
(without sign ambiguity) allow one to extract (f>s without any ambiguity'*. This then determines \Tto\. 



CDF and D0 have measured a*; directly and the average of their measurements is given by [33| 

all = -0.0115 ±0.0061 . 
If we take lAF^I = 0.075 ± 0.04, as given by CDF we obtain 

t&ncj), = -2.72 ±2.05 , iFfal = 0.11 ± 0.05 ps"* 

while lAF^I = O.lGStomt as given by D0 yields 



(51) 
(52) 



tan (A, = -1.25 ±0.83 



FJ2I = 0.13 ± 0.05 ps" 



(53) 



Although |Ff2l and 0s can significantly deviate from their SM predictions [Eq. ([5])], both of them are consistent with 
the SM within the error bar. Note that in the above numerical analysis we do not consider the negative solution for 
AFs, which introduces a sign ambiguity in the extraction of 4>s- It is clear that improved measurements of both a^; 
and AFs are essential in order to understand the underlying physics of B^-B^^ mixing and the width difference. 

In this paper we have focused on methods for measuring AFs using three-body decays. In practice, this will be 
carried out as follows. For deflnitiveness, consider the decays B'^{B'^) — >■ D^pK'^ K'^ . Generalising Eq. ((29)) to 
B^{B^) D'^pK^K^ , the time-dependent untagged differential decay distribution is given by 



untagged 



{D%pK+K-,t) 



(PT{B\ 







D%pK+K- 



(fTiB'i 



Dl;pK+K- 



ds^ds ds~^ds 
[A^f'" cosh(AFst/2) + Af,,^^ sinh(AFst/2)] , 



(54) 



where A^f^^ and Afi^^ are defined in Eq. ([5n| . Neglecting terms of order (AFs/F^)^ and higher, the time-integrated 
differential untagged decay distribution is given by 



dtTuntagged{D(jpK K ,t) — 



iDKK 
^ch 



2A 



DKK 
sh 



AFs 



(55) 



For a single resonance, say 



aDKK 

A DKK 

^sh 



Re 



-2i/3, 1^012 



2 —2i'y 

f e ^ 



(56) 



As discussed in the previous section, A^j^^ is fully known from the CP-averaged branching fraction of the intermediate 



resonance 4>. Once we have enough precision, a fit to the distribution given by Eq. (j54p or (j55|) allows one to obtain 
AFs and the various coefficients of ji^^P (which yields 2/3s). Such a fit will not allow the determination of the sign of 
AFs or cos(/)s, but Eq. (|5(I1) can still be used to obtain 0s (with a twofold ambiguity) and jF^jj. 

However, the above fit, though possible, is made difficult due to the requirement of having to simultaneously extract 
AFs and the components of A^j^^ . Given this, we would rather propose an alternative procedure. Referring again 
to Eq. (P^ . the time-dependent tagged differential decay distribution is given by 



- taggei 



d^T{B 







K- 



K- 



-r,t ^A^KK cos(Amst/2) - Af^^ sin(Amsi/2)] 



(57) 



where A^^^ and A^^^ are defined in Eqs. ((30)) . (|37)) and p4)) . From a fit to the above distribution, one can extract 
only the coefficients of different bilinears in A^^^ and A^^^ , since Atos is known. Thus, this fit straightforwardly 
gives information regarding A^^^ and A^^^ . As discussed in the previous section, from Af^^ alone we can extract 
2/?s and 7 + 2/3s ± S^, without any ambiguity, and 7 with the ambiguity [7, tt + 7]. This permits the reconstruction 
of [Eq. (HH)]. That is, all the coefficients of jF^p in Af/^ can be obtained from a fit to Eq. §7}. With this 



* Knowledge of tan<pB gives (ps with a twofold ambiguity, </>s •<-> tt + 0s. However, a"^ determines sin 0s, which allows one to differentiate 
(ps and 7r -t- 0s ■ 
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knowledge, there is only one unknown in Eq. ([54)) or (|55l) - AFg - and this can be determined by a fit. This may be 
a somewhat simpler procedure. Once we are able to measure AF^, then, along with a^j and AM^, Eq. ([50)1 can be 
used to obtain the CP phase (ps and iFfj]. 

The above analysis is also applicable to the decays B^{B'^) — > DfK^n^. However, as was discussed in Sec. IIII Cl 
for such decays all the trigonometric functions in A^j^ are not fully known - the only known functions are those 
appearing as the coefficients of \miFiF*) or Re(FiFj*) (i ^ j). Therefore, for such decays we can use Eq. (|54|) to fit 
AFg, but we need at least two interfering resonances, and only the terms proportional to lTtv{FiF*) or Re(FiFj*) are 

useful [see ^f^'^" in Eq. (p])]. 

V. CONCLUSIONS 

It is well known that the weak phase of B^^-B^ mixing is very small in the SM: 2/3^ ~ 0. If this quantity is measured 
to be significantly different from zero, this is a smoking-gun signal of new physics (NP). However, in general we would 
like more information from such a measurement. For instance, in order to distinguish among potential NP models, 
it is important to have an unambiguous determination of 2/3^. Similarly, although the width difference AF^ between 
the two Bs mass eigenstates is positive in the SM, it can take either sign in the presence of NP. Ideally, a method 
probing B^-B^ mixing which relies on a nonzero AFg should be able to remove its sign ambiguity. Finally, although 
it is usually assumed that NP contributes only to it has been shown that NP contributions to Ffj can also be 

important. In order to explore this possibility, it is necessary to measure the CP phase 0s and |F^2l- 

In this paper, we examine a variety of methods of measuring B^-B'^ mixing with an eye to addressing the above 
issues. We look at two- and three-body Bg decays with b cus and b ucs transitions, concentrating on those final 
states which are accessible to both i?° and 5° mesons (so that there is indirect CP violation). The time-dependent 
decay rates include both Amst and AFs^ terms. 

We begin with a review of B'^{B'^) — >■ DfK^ decays. Considering sizeable AFs, we find that this method allows 
the extraction of 2/3^ + 7 with a fourfold ambiguity. We then turn to B'^{B'^) — )• D^pcj) decays, where -D^'p ^ 
eigenstate. Here we find that 2/3^ and 27 can each be determined up to a twofold ambiguity. Here, the ambiguity is 
due to the unknown sign of AF^. Therefore, once we are able to resolve the sign ambiguity in AF^ by some other 
means, the B^{B^) — ^ D^p(f> decays are useful to measure 2/3s and 27 without any ambiguity. 

In order to resolve the sign ambiguity in AF^, and to reduce the discrete ambiguity in the measurement of 2/3s and 
7, it is necessary to turn to Dalitz-plot analyses of three-body decays. We begin with _B^(i?°) — > DfK^n^. We find 
that it is possible to obtain 2/3s -|- 7 with a twofold ambiguity, and to remove the sign ambiguity in AF^ (for this, it 
is not necessary to determine (ps). The most promising method involves the decays B'^(B'^) — > DqpKK, in which all 
issues can be resolved. We find that 2/3s can be obtained without any ambiguity, and at the same we can remove the 
sign ambiguity in AFc,. In addition, 7 can be determined up to a twofold ambiguity. 

Finally, all such decays allow the extraction of AFg directly from a fit to the time-dependent untagged differential 
decay rate distribution. Given the measurements of AM^, the semileptonic asymmetry a^;, and AF^, the CP phase 
(j)s and |Ff2l can be obtained. In the case of three-body decays the coefficients of sinh[AFst/2] and cosh[AFst/2] can 
be found, either fully or partially, from a fit to the time-dependent tagged differential decay rate distribution. (Of 
the several three-body decays that we discuss, the decays i?s(i?°) — D^pKK are the most promising, since in such 
decays these coefficients can be fully reconstructed from this fit.) Therefore, in three-body decays the only unknown 
in the untagged rate distribution is AF^. This makes the fit considerably simpler. 
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